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■ ABSTRACT 

<1> ' In this paper we study the role played both by novae and single stars in enriching 

the interstellar medium of the Galaxy with CNO group nuclei, in the framework of a 
£Nj ' detailed successful model for the chemical evolution of both the Galactic halo and disc. 

First, we consider only the nucleosynthesis from single low-mass, intermediate-mass 
and massive stars. In particular, the nucleosynthesis prescriptions in the framework of 
the adopted model are such that: i) low- and intermediate-mass stars are responsible 
for the production of most of the Galactic 12 C and 14 N; ii) massive stars produce 
the bulk of the Galactic 16 0; Hi) 13 C and 17 originate mostly in intermediate-mass 
stars, with only a minor contribution from low-mass and massive stars. In this context, 
we show that the behavior of the 12 C/ 13 C, 14 N/ 15 N and 16 0/ 17 isotopic ratios, as 
^ ■ inferred from observations, can be explained only allowing for a substantial revision 

| of the available stellar yields. On the other hand, the introduction of nova nucleosyn- 

f^i . thesis allows us to better explain the temporal evolution of the CNO isotopic ratios 

in the solar neighbourhood as well as their trends across the Galactic disc. Once all 
the nucleosynthesis sources of CNO elements are taken into account, we conclude 
that 13 C, 15 N and 17 are likely to have both a primary and a secondary origin, 
in contrast to previous beliefs. We show that, when adopting the most recent 17 
yields from intermediate-mass stars published in the literature so far, we still get a 
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too large solar abundance for this element, a problem already encountered in the past 
by other authors using different yield sets. Therefore, we conclude that in computing 
the 17 yields from intermediate-mass stars some considerable sink of 17 is probably 
neglected. The situation for 15 N is less clear than that for 13 C and 17 0, mainly due 
to contradictory observational findings. However, a stellar factory restoring 15 N on 
^ . quite long time-scales seems to be needed in order to reproduce the observed positive 

gradient of 14 N/ 15 N across the disc, and novae are, at present, the best candidates 
for this factory. Given the uncertainties still present in the computation of theoreti- 
cal stellar yields, our results can be used to put constraints on stellar evolution and 
nucleosynthesis models. 

Key words: Galaxy: abundances - Galaxy: evolution - novae, cataclysmic variables 
- nuclear reactions, nucleosynthesis, abundances 



1 INTRODUCTION 

In the last 30 years, many studies have been devoted to the 
evolution of the CNO isotopic ratios (e.g., Audouze, Lequeux 
& Vigroux 1975; Vigroux, Audouze & Lequeux 1976; Dear- 
born, Tinsley & Schramm 1978; Tosi 1982; D'Antona & Mat- 
teucci 1991; Matteucci & D'Antona 1991; Prantzos, Aubert 
& Audouze 1996; see also Tosi 2000 for a recent reappraisal 
of the problem) . Isotopic ratios are generally not affected by 
physicochemical fractionation effects (but see, e.g., Sheffer, 



Lambert & Federman 2002). Therefore, they reflect rather 
faithfully the relevant production processes which occurred 
in stars of different masses and lifetimes. However, despite 
the considerable progress in the theories of stellar evolu- 
tion and nucleosynthesis, important questions related to the 
yields and the production sites of some of the CNO isotopes 
still remain open. In this context, chemical evolution models 
can be regarded as a powerful tool in order to discriminate 
among different sets of stellar yields and even different stellar 
factories. In fact, as discussed by Tosi (1988), model predic- 
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tions on elemental and isotopic ratios depend primarily on 
the adopted stellar nucleosynthesis, stellar lifetimes and ini- 
tial mass function (IMF) , rather than on the other uncertain 
parameters of the galactic evolution, like star formation and 
infall rates. 

Among the CNO group nuclei, 16 O is the best under- 
stood. It is a primary element, i.e., it is always synthesized 
starting from H and He in the parent star. When only the 
most reliable oxygen measurements are considered, a good 
agreement is found between theoretical predictions on 16 O 
evolution and observations (e.g., Chiappini, Romano & Mat- 
teucci 2003), which confirms that 16 O production is well un- 
derstood both qualitatively and quantitatively. The bulk of 
O comes from massive stars (m > 10 Mq). 

The production of 12 C is instead more uncertain: we 
know that it is synthesized as a primary element, but the 
exact amount restored into the interstellar medium (ISM) by 
stars of different masses is still uncertain. The 'old' view that 
12 C is mainly produced by low- and intermediate-mass stars 
(LIMS) has been challenged by several papers (Prantzos, 
Vangioni-Flam & Chauveau 1994; Carigi 2000; Henry, Ed- 
munds & Koppen 2000) suggesting that mass-loosing mas- 
sive stars could be the main source of carbon. The conclusion 
of those authors rests on the adoption of the metallicity de- 
pendent yields for massive stars from Maeder (1992), which 
are computed in presence of mass loss. On the other hand, 
very recently it has been demonstrated that, in order to re- 
produce simultaneously the features of the C/O vs. O/H di- 
agrams observed for galaxies of different morphological type, 
as well as the [C/Fe] vs. [Fe/H] and [C/O] vs. [Fe/H] dia- 
grams for dwarf stars in the solar vicinity, the bulk of 12 C 
should originate from LIMS (Chiappini et al. 2003). Indeed, 
the carbon yields due to mass loss have been revised down- 
ward (Meynet & Maeder 2002a). 

The nature of 14 N is still controversial: it should be a 
typical secondary element, but with a primary component as 
well. The secondary component is mainly produced through 
the CN cycle in stars of all masses, thus its amount de- 
pends on the 12 C present in the star at the time of its birth. 
Secondary nitrogen can also be produced in the ON cycle 
by transformation of 16 O, but at a much slower rate. The 
primary component is instead produced starting from fresh 
carbon generated by the parent star, provided there is some 
mixing mechanism which transports this newly synthesized 
carbon in a hydrogen-burning region, where the CNO cy- 
cle can convert it to nitrogen. In intermediate-mass stars 
(2.5-3.5 < tu/Mq < 6-8), primary 14 N is produced when 
hot bottom burning (HBB) takes place at the base of the 
star convective envelope, during the asymptotic giant branch 
(AGB) phase (Renzini & Voli 1981). The amount of primary 
nitrogen depends mostly on the adopted value of the mixing 
length parameter determining the depth of the convective 
stellar envelope, which is very uncertain (see a discussion in 
van den Hoek & Groenewegen 1997). Primary 14 N produc- 
tion might take place also in massive stars (e.g., Matteucci 
1986; Meynet & Maeder 2002b). According to Meynet & 
Maeder (2002b), the mechanism which allows the freshly 
made carbon to be turned into primary nitrogen is stellar 
axial rotation, and primary nitrogen production in massive 
stars is favored at low metallicity. In fact, in this case (ro- 
tating stellar models at low Z), the stars loose less angular 
momentum, thus they rotate faster. Simultaneously, they 



are more compact, so that differential rotation and shear 
mixing arc stronger. Moreover, the H-burning shell has a 
much higher temperature and is thus closer to the core, 
which favours mixing between the two. However, the pri- 
mary nitrogen yields obtained by Meynet & Maeder from 
their low-metallicity massive star models do not contribute 
to a significant chemical enrichment at very low metallicities 
(see Chiappini et al. 2003). 

As far as the minor isotopes are concerned, it has been 
claimed that 13 C and possibly also 17 are mainly formed in 
the external regions of stars in the red giant branch (RGB), 
planetary nebula (PN) and supernova (SN) phase (Audouze 
et al. 1975; Dearborn et al. 1978). Marigo, Bressan & Chiosi 
(1996) explored the formation of a pocket of primary 13 C 
during the AGB phase of intermediate-mass stars. 

The thermonuclear runaway (TNR) responsible for the 
nova outbursts has also been identified as a promising chan- 
nel for the synthesis of 13 C, 15 N, and 17 (Starrfield et al. 
1972; Starrfield, Sparks & Truran 1974), eventually being 
responsible for most of the Galactic abundances of 13 C and 
17 (Jose & Hernanz 1998; Woosley et al. 1997). Indeed, 
classical novae, binary systems consisting of a carbon-oxygen 
or oxygen-neon-magnesium white dwarf accreting hydrogen- 
rich matter from a main-sequence companion, sporadically 
inject nuclearly processed material into the ISM. The TNR, 
responsible for the explosion causing the ejection of almost 
the whole previously accreted envelope, leads to the synthe- 
sis of some rare nuclei during the nova outburst. The most 
recent nova hydrodynamical models confirm large overpro- 
duction factors with respect to the solar abundances for 7 Li, 
13 C, 15 N and 17 (Jose & Hernanz 1998). Therefore, al- 
though novae are processing only a minor fraction of the 
total matter in a galaxy, nevertheless they could produce 
important amounts of the species with the highest overpro- 
duction factors. In the past, it has been suggested that the 
bulk of 15 N production should come from novae (Audouze et 
al. 1975; Dearborn et al. 1978; Matteucci & D'Antona 1991); 
Type II SNe might also contribute, although to a lower 
extent (Audouze et al. 1977; Timmes, Woosley & Weaver 
1995). Recent detections of extragalactic 15 N seem to sup- 
port the idea that 15 N production is due to rotationally in- 
duced mixing of protons into the helium-burning shells of 
massive stars (Chin et al. 1999 and references therein). A 
small contribution from Type la SNe could also be possible 
(Clayton et al. 1997; although see Nomoto, Thielemann & 
Yokoi 1984). 

The aim of this paper is to present some theoretical 
results concerning the evolution of the 12 C/ 13 C, 14 N/ 15 N, 
and 16 0/ lr O isotopic ratios in the solar neighbourhood as 
well as along the Galactic disc by taking advantage of a 
complete model for the chemical evolution of the Galaxy in 
which the nucleosynthesis from stars of different masses is 
included in a detailed way. In particular, we include in the 
model the nucleosynthesis from novae and try to establish 
whether they can actually be regarded as important contrib- 
utors to the Galactic abundances of the CNO isotopes 13 C, 
15 N, and ir O. For each of these stable nuclei, we make an 
attempt in identifying the main contributors to the chemi- 
cal enrichment by means of a comparison with the available 
observational data. 

The paper is organized as follows. In Section 2 we re- 
port on the status of the observations. In Section 3 we de- 
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scribe the adopted chemical evolution model. In Section 4 
we present model results. Finally, our conclusions are drawn 
in Section 5. 



2 OBSERVATIONS 

2.1 Galactic 12 C/ 13 C abundance ratio 

The Solar System 12 C/ 13 C ratio is ( 12 C/ 13 C) = 89 ± 2 
(Cameron 1982), higher than the estimates for the local ISM: 
( 12 C/ 13 C)ism = 47.3± |-f (from CN; Crane & Hegyi 1988); 
77 ± 3 (from CH+; Stahl et al. 1989); 62 ± 4 (from CO; 
Langer & Penzias 1993); 77 ± 7 (average from H2CO and CO 
data for sources near the Sun; Wilson & Rood 1994); 58 ± 12 
(from atomic C; Keene et al. 1998); 58 ± 6 (from atomic C + ; 
Boreiko & Betz 1996); 69 ± 15 (from the solid state CO2; 
Boogert et al. 2000)*. This clearly indicates a decrease in the 
last 4.5 Gyr, which is exactly what one expects for a primary 
to secondary elemental ratio (according to Talbot & Arnett 
1974, a chemical species is defined as primary or secondary 
according to whether the yields, i.e., the mass fractions of 
a star ejected in the form of the newly produced element, 
are insensitive, or sensitive, to the original metallicity of 
the star; as a consequence, the abundance of a secondary 
element increases more steeply in time than that of its seed 
nuclei) . 

The Galactic 12 C/ 13 C gradient gives a consistent pic- 
ture: A( 12 C/ 13 C)/AR G = (4.5 ± 2.2) kpc" 1 (weighted fit 
to the 12 C0 2 / 13 C0 2 ratio as a function of Galactocentric 
radius - Boogert et al. 2000). The gradient found from the 
solid state agrees well with gas phase studies (Tosi 1982; 
Wilson & Rood 1994). However, solid state determinations 
are more reliable than gas phase studies, because the col- 
umn densities can be determined without uncertain radia- 
tive transfer effects. A limitation of the Boogert et al. anal- 
ysis is the paucity of solid CO2 observations at low Galactic 
radii (3-6 kpc). Moreover, the considerable scatter in the 
12 C/ 13 C ratio at radii larger than 6 kpc prevents from identi- 
fying a clear trend in the outer disc. However, the decrease of 
the ratio towards the Galactic centre is apparent and can be 
taken as a clear signature of an important secondary com- 
ponent to 13 C nucleosynthesis. This decrease is confirmed 
also by 12 C/ 13 C ratios derived from 12 CN, 13 CN observa- 
tions, which however tend to be lower than those derived 
from other molecules (Savage, Apponi & Ziurys 2001). 

2.2 Galactic 14 N/ 15 N abundance ratio 

The 14 N/ 15 N ratio in our Galaxy has been studied by Dah- 
men, Wilson & Matteucci (1995), who analysed high qual- 
ity data of HCN double isotopomer ratios for 11 warm, 
dense molecular clouds in the Galactic disc. After fitting 
12 C/ 13 C ratios from published data as a function of Galac- 
tocentric distance, the 14 N/ 15 N ratios were obtained. A lin- 
ear regression using these ratios and previously published 
data allowed the determination of the Galactic gradient: 
A( 14 N/ 15 N)/AR G = (19.7 ± 8.9) kpc" 1 , shallower than 
that found by Tosi (1982). 

* A value of ( 12 C/ 13 C)ism between 60 and 80 is the preferred 
one. 



The Solar System 14 N/ 15 N ratio is ( 14 N/ 15 N) Q = 270 
(Anders & Grevesse 1989). For the local ISM, values as high 
as ( 14 N/ 15 N)ism = 450 ± 22 (Wilson & Rood 1994) or as low 
as ( 14 N/ 15 N)i SM = 237^21 (Lucas & Liszt 1998) have been 
reported. The former is in contradiction with the existence 
of a positive gradient across the Galactic disc since, from a 
theoretical point of view, if 14 N/ 15 N increases in time and 
hence with metallicity, it should also increase with decreas- 
ing Galactocentric distance, owing to the higher metallicity 
content attained in the innermost Galactic regions. 

2.3 Galactic 16 0/ 17 abundance ratio 

The 16 0/ 17 abundance ratio has decreased from its solar 
value of ( 16 O/ 17 O) = 2465 (Anders & Grevesse 1989) to 
its present local value of ( 16 0/ 17 0)ism = 1900 ± 200, an 
unmistakable sign of secondary production for 17 0. 

The secondary origin of 17 O is confirmed by the behav- 
ior of the le O/ 17 gradient along the Galactic disc, which 
decreases with decreasing distance from the Galactic centre 
( le O/ 17 ratios along the Galactic disc are obtained by mul- 
tiplying the corresponding 16 0/ 18 ratios by the ISM value 
of 18 0/ 17 = 3.5 - see Prantzos et al. 1996 and references 
therein). 



3 THE CHEMICAL EVOLUTION MODEL 
3.1 Basic aspects 

The two-infall model for the chemical evolution of the 
Galaxy (Chiappini, Matteucci & Gratton 1997; Chiappini, 
Matteucci & Romano 2001) is adopted, with the aim of re- 
assessing the problem of whether novae significantly con- 
tribute to the Galactic chemical evolution of the minor CNO 
isotopes. 

In the framework of the adopted model, the Galaxy 
forms out of two main infall episodes. During the first one, 
the bulge, inner halo and thick disc are built up, on a very 
short time-scale (~ 1 Gyr). During the second one, the thin 
disc is formed, on a much longer time-scale (being ~ 7 Gyr at 
the solar ring and increasing with increasing Galactocentric 
distance) out of matter of primordial chemical composition 
plus traces of halo gas. 

For a complete description of the model basic assump- 
tions and equations we address the interested reader to Chi- 
appini et al. (1997, 2001, 2003). Here we only briefly out- 
line how novae have been included in the model (see also 
DAntona & Matteucci 1991 and Romano et al. 1999). The 
rate of formation of nova systems at a given time t is com- 
puted as a fraction a of the rate of formation of white dwarfs 
(WDs) at a previous time t — At: 

,8 

-Rnovao(t) = a / tp(t -T m - At)ip(m)dm, (1) 
Jo.s 

where r m is the lifetime of the star of mass m, ip(t) is the 
star formation rate and ip(m) is the IMF. The Scalo (1986) 
IMF is assumed. All the stars between m = 0.8 Mq and m 
— 8 Mq end up as WDs. At is a suitable delay-time which 
guarantees the cooling of the WD at a level that ensures 
a strong enough nova outburst. Here we assume At = 2 
Gyr, which is a typical value for novae (see also Romano 
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et al. 1999 and D'Antona & Mazzitelli 1982). The parame- 
ter a represents the fraction of WDs which belong to nova 
systems and its value (constant in time) is fixed by the re- 
quest of reproducing the current rate of nova outbursts in 
the Galaxy. Unfortunately, this is a free parameter. Here we 
set a ~ 0.01. The rate of nova outbursts is computed from 
the rate of formation of nova systems by assuming that each 
nova suffers roughly 10 4 outbursts during its life (Bath & 
Shaviv 1978). A value of a between ~ 0.01 and ~ 0.02 leads 
to -Routbursts(tGai) ~ 20 - 30 yr _1 , which has to be compared 
with the value inferred from scalings from extragalactic nova 
surveys, iO,ursts(<Gai) = 15-50 yr" 1 (Delia Valle & Livio 
1994; see also Shatter 1997). 

Since the evolution of the CNO isotopic ratios (as well 
as that of the other elemental ratios) depends mainly on the 
nucleosynthesis prescriptions, rather than on the parame- 
ters of the galactic evolution like star formation and infall 
rates, chemical evolution models can put constraints mainly 
on the nucleosynthesis and time-scales for element produc- 
tion. The validity of this statement rests on the fact that 
any variation in the star formation and/or infall laws affects 
the main and the minor isotopes at the same extent. In the 
framework of a model where nova nucleosynthesis is taken 
into account, the situation is more complicated. In fact, in 
this case the evolution of the CNO isotopic ratios depends 
also on the parameters which enter the computation of the 
theoretical nova outburst rate. Since novae produce large 
amounts of the minor isotopes, whereas their contribution 
to the main ones is almost negligible, each variation in the 
parameters regulating the nova rate does not cancel out in 
the CNO isotopic ratios. One should always keep in mind 
this when discussing predictions from models in which nova 
nucleosynthesis is included. We will come back on this issue 
in Section 4.3. 

3.2 Nucleosynthesis prescriptions 

The nucleosynthesis prescriptions adopted here are from 
i) van den Hoek & Groenewegen (1997) and Ventura, 
D'Antona & Mazzitelli (2002) for LIMS; ii) either Woosley 
& Weaver (1995) or Nomoto et al. (1997) for Type II SNe; 
Hi) Thielemann et al. (1993) for Type la SNe; iv) Jose & 
Hernanz (1998) for novae. 

It is worth emphasizing that metallicity dependent 
yields are available from a limited number of studies, espe- 
cially in the range of massive stars (m > 10 Mq). Moreover, 
nucleosynthesis studies are usually restricted to specific mass 
ranges (e.g., m ~ 0.9-8.0 Mq - van den Hoek & Groenewe- 
gen 1997; m ~ 2.5-6.0 Mq - Ventura et al. 2002; m ~ 
11-40 M e - Woosley & Weaver 1995; m ~ 13-70 Mq - 
Nomoto et al. 1997) and/or do not deal with some specific 
chemical species, so that it is neither possible to homoge- 
neously cover the mass spectrum over which stars distribute 
(m ~ 0.1-100 Mq) nor to treat all the relevant chemical 
species by adopting a single stellar nucleosynthesis study. 

For instance, van den Hoek & Groenewegen (1997) give 
metallicity dependent yields of 4 He, 12 C, 13 C, 14 N and 16 
for stars in the mass range m ~ 0.9 - 8.0 Mq , but do not pro- 
vide the yields of 17 O. Therefore, one must complete their 
grid of stellar yields by means of 17 O yields coming from 
some other study. We adopt the lr O yields recently com- 
puted by Ventura et al. (2002) for stars in the mass range 



m ~ 2.5-6.0 Mq, for different initial chemical compositions. 
Since their yields do not extend to metallicities higher than 
Z = 0.01, we must extrapolate them to supersolar metallic- 
ities. This makes the modelisation of 17 evolution in the 
inner disk inaccurate. Stars with m < 2.0 Mq are not 17 
producers, since they do not go through the NO cycle during 
which this element is synthesized. Therefore, we set to zero 
the lr O yields from stars in the mass range m ~ 1.0-2.0 
Mq. We also set to zero the 15 N production from LIMS. In 
fact, the net yield of this element from LIMS is negative, 
due to the effect of the first and second dredge-up for stars 
in the mass range 0.8-3.5 Mq and to HBB for higher mass 
stars (Marigo 2001). 

Among the input parameters of nova hydrodynamical 
models with a deep influence on the nova nucleosynthesis 
there is the chemical composition of the H-rich envelope ac- 
creted by the WD from the main-sequence companion which 
fills its Roche lobe. The problem of the chemical composition 
of nova envelopes is complex and far from being understood. 
The matter transferred from the companion is assumed to 
be solarlike and is mixed in a given fraction with the out- 
ermost shells of the underlying core (e.g., Politano et al. 
1995; Jose & Hernanz 1998). This is done in order to get 
the enhanced CNO or ONeMg abundances required both to 
power the explosion and to account for the spectroscopic 
abundance determination. It is not clear how the composi- 
tion of the nova ejecta changes during the evolution of the 
Galaxy, as a function, e.g., of the chemical composition of 
the matter accreted from the WD companion. Dearborn et 
al. (1978) assumed that the amount of 5 N ejected per nova 
outburst scales with the initial 12 C + 13 C + 14 N abundance 
of the star. This means that 15 N from novae is an element 
of purely secondary origin. On the contrary, la C and 15 N 
from novae were assumed to be of purely primary origin by 
D'Antona & Matteucci (1991) and Matteucci & D'Antona 
(1991). However, as noticed by these authors, a primary ele- 
ment mainly produced in long-lived progenitors still behaves 
like a secondary one from the point of view of galactic chem- 
ical evolution, owing to the delay with which it is restored to 
the ISM (see also Matteucci & Francois 1989). In Section 4.3 
we will show results obtained by assuming two limiting cases 
for 13 C, 15 N and 17 production during nova outbursts: in 
the first case, the yields of 13 C, 15 N and lr O from novae are 
assumed to be always the same, at each time of Galaxy's 
evolution; in the second case, they are scaled with the abun- 
dance of the seed nuclei in the star. 

The nucleosynthesis from novae is taken into account in 
the chemical evolution code as follows. We consider that 30% 
of novae arise from systems containing ONe WDs and 70% 
from systems containing CO WDs (e.g., Gehrz et al. 1998) 
and compute the abundance (in mass fraction) of element i 
in the ejecta of the typical nova as: 

(Xi) = 0.3 x (Xi) 

ONc nova + 0.7 x {Xi) 

C O nova j (2) 

where (Vi)oNonova and (A;)conova are the average abun- 
dances of the element i in the ejecta of ONe novae and CO 
novae, respectively (Romano et al. 1999). These quantities 
are obtained in turn by averaging the results of two grids 
of 7 evolutionary sequences each computed by Jose & Her- 
nanz (1998) for ONe and CO WDs, respectively (see their 
tables 3 and 4). Once the elemental abundances in the ejecta 
of the typical nova have been estimated, we still need to as- 
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Table 1. Nucleosynthesis prescriptions. 



Model 


clement 


LIMS 


massive stars a 


novae: ejecta^ 


novae: s/p 


1 


12 C 
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- 




13 C 


vdHG97 (r, const ) 


N97 
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- 
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TS.TQ7 
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no 
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VDM02/6 


N97 
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2 


12 C 


vdHG97 (r, const ) 


N97 


6.7 x 10- 3 M Q 


P 




13 C 


no 


no 


1.2 x 10~ 2 Mq 


P 




14 N 


vdHG97 (r, const ) 


N97 


1.7 X 10~ 2 Mq 


P 




15 N 


no 


no 


1.4 x ltr 3 Mq 


P 




16o 
w 


™nby / V^const ) 


]\TQ7 




P 




17 o 


no 


no 


1.4 X 10^3 Mq 


P 


2s 


12 C 


vdHG97 (rfconst) 


N97 


6.7 x 10- 3 M 


P 




13 C 


no 


no 


1.5 X 1(T 2 Mq 


s 




14 N 


vdHG97 (ryconst) 


N97 


1.7 X 1(T 2 Mq 


P 




15 N 


no 


no 


3.8 x 10~ 3 Mq 


s 




16 o 


van^y i v '/const ; 


NQ7 
in y / 




p 




17 o 


no 


no 


1.8 X 10^3 Mq 


s 


3 


12 C 


vdHG97 (fjconst) 


N97 


6.7 x 10^3 Mq 


p 




13 C 


vdHG97 ()jc„„ s t) 


N97 


5.3 x 10^3 Mq 


p 




14 N 


vdHG97 ()jc„„ s t) 


N97 


1.7 x 10~ 2 M 


p 




15 N 


no 


no 


1.4 X 10^3 Mq 


p 




16q 


vunuy/ ^ '/const ) 


N97 


9^*1 fl^ 1 Mrs 


p 




17 o 


VDM02/10.5 


N97 


6.4 x 10- 4 M 


p 


3s 


12 C 


vdHG97 ()jc„„ s t) 


N97 


6.7 x 10-3 Mq 


p 




13 C 


vdHG97 (r; const ) 


N97 


1.5 x 10- 2 Mq 


s 




14 N 


vdHG97 ()jc„ M t) 


N97 


1.7 x 10- 2 Mq 


p 




15 N 


no 


no 


3.8 x 10-3 Mq 


s 




16 o 


vunuy/ ^ '/const ; 


N97 




p 




17 o 


VDM02/10.5 


N97 


8.8 x 10- 4 Mq 


s 


3n 


12 C 


vdHG97 (?? v ar) 


WW95 


6.7 x 10-3 Mq 


p 




13 C 


vdHG97 (r/var) 


WW95 


1.5 x 10- 2 M 


s 




14 N 


vdHG97 (»/var) 


WW95 


1.7 x 10- 2 Mq 


p 




1B N 


no 


no 


3.8 x 10-3 Mq 


s 




16 r\ 


vdHG97 (»/var) 


WW95 


2.5 X 10 Mq 


p 




17 o 


VDM02/10.5 


WW95 


8.8 x 10- 4 Mq 


s 


3m 


12 C 


vdHG97 (J7 var ) 


WW95 


6.7 x 10-3 Mq 


p 




13 c 


vdHG97 (»7var) 


WW95 


1.5 x 10- 2 Mq 


s 




14 N 


vdHG97 (»7var) 


WW95 


1.7 x 10- 2 Mq 


p 




15 N 


no 


no 


3.8 x 10-3 M 


s 




16 Q 


vdHG97 (»7var) 


WW95 


2.5 x 10- 1 Mq 


p 




17 o 


VDM02/15 


WW95 


8.8 x 10- 4 Mq 


s 


a The yield 


s of 12 C have been multiplied by a factor of 3 


in the range 40- 


100 Mq (see 



arguments in Chiappini et al. 2003). 

b The mean masses ejected in form of 12 C, 13 C, 14 N, 15 N, ie O and 17 by each nova 
during its overall lifetime arc: 6.7 X lO -3 Mq, 1.5 X lO -2 Mq, 1.7 X 10 -2 Mq, 5.2 X 
IO- 3 Mq, 2.5 X lO -1 Mq and 3.2 X lO -3 Mq, respectively. These quantities are obtained 
by averaging the yields over the 14 evolutionary sequences given by Jose & Hernanz (1998) 
and by assuming that the mean mass ejected in a single burst is ~ 2 X 10~ 5 Mq. They refer 
to Z = Zq, and arc computed by taking into account that 30% of nova systems contain 
ONe WDs, while the remaining contain CO WDs. The quantities listed in column 5 are 
those actually needed in order to reproduce the solar abundances of 13 C, 15 N and 17 as 
given by Anders & Grevesse (1989). It is seen that in many cases it is necessary to lower 
the quantities given by the average described above. 
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Figure 1. a) Temporal evolution of the carbon isotopic ratio in the solar neighbourhood in the case of 13 C production from single stars 
alone (Model 1, dashed line), novae alone (Model 2, dots) or novae plus single stars (Model 3, continuous line). The model 12 C/ 13 C 
ratios 4.5 Gyr ago are normalised to their solar values. The solar abundances by mass of 12 C and 13 C predicted by the models are given 
at the bottom of the figure. The big dot is the observed meteoritic ratio; the vertical bar is the ratio observed in the local ISM (references 
can be found in the text). The minimum in the curves at around 1 Gyr is due to the transition between halo-thick disc and thin disc in 
the two-infall model, b) Theoretical present-day carbon isotopic ratio across the Galactic disc in the case of 13 C production from single 
stars alone (Model 1, dashed line), novae alone (Model 2, dots) or novae plus single stars (Model 3, continuous line) vs. observations. 



sume a fiducial value for the mass ejected during each nova 
outburst. This is a highly uncertain quantity. The observa- 
tionally derived mass distribution of nova ejecta indicates 
that the data peak close to a few 1O~ 4 M (Delia Valle 2000 
and references therein), but this estimate could be seriously 
challenged if the filling factor of the ejected shells is closer 
to 10~ 2 - 10" 5 rather than ~ 0.1- 1 (e.g., Shara et al. 1997). 
We assume M cj ~ 2 x 10~ 5 M®, in agreement with the 
predictions of the mass ejected during a nova outburst from 
nova hydrodynamical models, and compute the mass of the 
element i ejected during the whole nova evolution as: 

M? = (Xi) x M ci x n, (3) 

where n = 10 4 is the number of nova outbursts suffered on 
an average by each nova system. 

In columns 3 and 4 of Table 1 we list the nucleosynthesis 
prescriptions for 12 C, 13 C, 14 N, 15 N, 16 and 17 from LIMS 
and massive stars adopted in our different models. Models 
are labelled as given in column 1. The average masses of 12 C, 
13 C, 14 N, 15 N, 16 and 17 ejected by each nova system 
during its overall evolution are listed in column 5. In some 
cases, they differ from the average values obtained with the 
procedure described above (Eq. 3), since these latter tend 
to overestimate the solar abundances of 13 C, 15 N and ir O, 
while in this table we want to list the contributions we ac- 
tually used in order to reproduce the observations. For each 
model, column 6 indicates if the CNO production from novae 
is treated as primary or secondary. When the CNO produc- 
tion is treated as secondary, the quantities listed in Table 1, 
corresponding to accreted envelopes of solar chemical com- 



position, have been scaled according to the abundances of 
the seed nuclei predicted by the model at each time. 



4 RESULTS 

4.1 13 C production from LIMS 

It is well known that LIMS produce 13 C, both as a pri- 
mary and secondary element. 13 C is also produced by mas- 
sive stars, although in smaller amounts. In Fig. 1 we show 
as a dashed line the results of a model (Model 1, Table 1) 
in which 13 C is made only by single stars. Nucleosynthesis 
prescriptions for LIMS and massive stars are from van den 
Hoek & Groenewegen (1997 - their standard case in which 
the mass loss parameter for stars on the AGB, t^agb , re- 
mains constant and equal to 4) and Nomoto et al. (1997), 
respectively. It is found that although single stars alone can 
explain the solar abundance of 13 C (Table 2), they fail in 
accounting for the observed decrease of the 12 C/ 13 C ratio 
in the solar neighbourhood in the last 4.5 Gyr (Fig. 1, left 
panel). Conversely, they predict a roughly constant 12 C/ 13 C 
ratio in the past ~ 12 Gyr. The same conclusions still hold 
if we adopt the van den Hoek & Groenewegen (1997) yields 
computed with ijagb = 1 for Z — 0.001 and ijagb = 2 for 
Z = 0.004 for LIMS, and the Woosley & Weaver (1995) 
yields rather than the Nomoto et al. (1997) ones for massive 
stars. This is due to the fact that in all these models the 
bulk of 13 C comes from intermediate-mass stars suffering 
HBB (m > 4 M Q , r m < 180 Myr). Therefore, changing the 
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Table 2. Predicted and observed abundances by mass at the time of Sun formation (£q = 9.5 
Gyr). 

12 C 13 C 14 N 15 N 16 17 Q 

Model 1 

1.47 x 1CT 3 2.21 x 1CT 5 6.92 x 10~ 3 3.87 x 10~ 6 

Model 2 

1.49 x 1CT 3 3.87 x 1CT 5 0.96 x 10~ 3 3.87 x 10~ 6 7.72 x 10~ 3 3.89 x 10~ 6 

Model 2s 

1.49 x 1CT 3 1.76 x 1CT 5 0.95 x 10~ 3 3.87 x 10~ 6 7.71 x 10~ 3 3.85 x 10~ 6 

Model 3 

1.49 x 10~ 3 3.97 x 10~ 5 0.95 x 10~ 3 3.88 x 10~ 6 7.71 x 10~ 3 4.00 x 10~ 6 

Model 3s 

1.49 x 10~ 3 3.97 x 10~ 5 0.95 x 10~ 3 3.87 x 10~ 6 7.71 x 10~ 3 3.99 x 10~ 6 

Model 3n 

1.86 x 10~ 3 4.67 x 10~ 5 1.00 x 10~ 3 4.99 x 10~ 6 6.12 x 10~ 3 4.41 x 10~ 6 

Model 3m 

1.86 x 10~ 3 4.67 x 10~ 5 1.00 x 10~ 3 4.99 x 10~ 6 6.12 x 10~ 3 4.44 x 10~ 6 

Observed — Anders & Grevesse (1989) 

3.03 x 10~ 3 3.65 x 10~ 5 1.11 x 10" 3 4.36 x 10~ 6 9.59 x 10~ 3 3.89 x 10~ 6 

Observed - Holwcgcr (2001) 
3.31 x 10~ 3 0.85 x 10~ 3 6.24 x 10~ 3 

Observed - Allende Prieto et al. (2001, 2002) 
2.09 x 10~ 3 5.56 x 10~ 3 



yields either in massive stars or in the low-metallicity range 
does not affect our main conclusions. 

During HBB, 13 C is mostly produced as a primary el- 
ement. However, a secondary component is present as well, 
depending on how much of the dredged-up 12 C is newly pro- 
duced inside the star and how much was already present at 
the time of the star birth (see van den Hoek & Groenewegen 
1997). The constancy of the 12 C/ 13 C ratio we predict in the 
last ~ 12 Gyr is mainly due to the fact that most of the 13 C 
is of primary origin, and to the fact that in the framework 
of our model 13 C production proceeds almost in step with 
that of 12 C. In order to solve the discrepancy with the obser- 
vations, one may argue whether the 13 C yields from single 
low-mass stars, which produce 13 C always as a secondary 
element, are underestimated. 

Indeed, in the last few years a great deal of work has 
appeared in the literature on whether standard stellar nu- 
cleosynthesis models are underestimating the 13 C content in 
the material ejected from LIMS. It has often been suggested 
that the progenitors of the PNe which are observed to have 
carbon isotopic ratios below the values expected from stan- 
dard AGB models must have undergone some non-standard 
mixing process, perhaps induced by rotation, during their 
red giant phase and/or AGB phase (Bachiller et al. 1997; 
Palla et al. 2000; Balser, McMullin & Wilson 2002). The 
existence of this non-standard mixing mechanism is also in- 
voked in order to reconcile predictions from standard Galac- 
tic evolutionary models with observations of 3 He (e.g., Galli 



et al. 1997; Charbonnel & do Nascimento 1998; Chiappini, 
Renda & Matteucci 2002). 

However, the extent of the 13 C enhancement as a func- 
tion of mass with varying metallicity is, to date, rather un- 
certain. In fact, the requirement of sampling low-mass ob- 
jects has often led to concentrating the observational efforts 
on metal-poor clusters. Charbonnel, Brown & Wallerstein 
(1998) found that field giants belonging to the old Galactic 
disc population ([Fe/H] ~ -0.6) can reach 12 C/ 13 C ~ 7. For 
more metal-rich stars, a reasonable estimate might be ~ 10 - 
15 (from observations of the carbon isotopic ratio in stars as- 
cending the RGB in the open cluster M 67; Gilroy & Brown 
1991). Measurements of 12 C and 13 C abundances in eleven 
bright giant members of the globular cluster u> Centauri 
spanning the metallicity range —1.8 < [Fe/H] < —0.8 do not 
reveal any real trend of the efficacy of the first giant branch 
mixing with metallicity (Smith, Terndrup & Suntzeff 2002). 
The mean value for this sample is ( 12 C/ 13 C) = 4.3±0.4, with 
nine stars sharing, within the errors, the equilibrium ratio 
of 3.5. Standard stellar evolution theory predicts a decrease 
of 12 C/ 13 C from ~ 90 to ~ 25 during the first dredge-up 
on the RGB for LIMS (Iben 1964; Iben & Renzini 1984). 
Therefore, an extra-mixing which further processes material 
during the RGB seems to be at work. 

The 12 C/ 13 C ratio in Galactic carbon stars has been 
observed by Lambert et al. (1986), Onhaka & Tsuji (1996, 
1999), Abia & Isern (1997) and Schoier & Olofsson (2000). In 
particular, Onhaka & Tsuji (1996, 1999) found 12 C/ 13 C ra- 
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Figure 2. 12 C/ 13 C vs. progenitor mass in the total ejecta of 
LIMS. Stellar models are both from van den Hoek & Groenewegen 
(1997) (upper panel) and Ventura et al. (2002) (lower panel), and 
have been computed for different values of the initial metallicity 
of the stars. The box represents the region occupied by the PNe 
detected in CO by Palla et al. (2000). Ventura et al. (2002) predict 
sistematically lower ratios than van den Hoek & Groenewegen 
(1997) in the mass range 3-6 Mq. However, their computations 
do not extend to the interesting 1-3 Mq range, where model 
predictions could be usefully compared to observations. 



tios about a factor of 2 to 3 smaller than those of Lambert et 
al. (1986) for the same starsT, whereas Abia & Isern (1997) 
found something falling in between the results obtained by 
Lambert et al. (1986) and Onhaka & Tsuji (1996). To solve 
the controversy, Schoier & Olofsson (2000) performed inde- 
pendent estimates of the 12 C/ 13 C ratio, using CO radio line 
emission from the circumstellar envelopes, for a sample of 
carbon stars showing large discrepancies between the sets 
of photospheric estimates. They used a non-LTE radiative 
transfer code. The 12 C/ 13 C ratios range from 2.5 to 90 and 
agree with those estimated in the photospheres by Lambert 
et al. (1986). Finally, in a recent paper, Abia et al. (2001) 
suggest the activation of slow circulation mechanisms below 
the convective border of the envelope for most of the ob- 
served stars, in order to explain the low observed 12 C/ 13 C 
values. In conclusion, it appears that at least a fraction of 
PN progenitors do have reduced 12 C/ 13 C ratios compared 
to standard model expectations, and that some extra-mixing 
should be required. 

In Fig. 2 we show the 12 C/ 13 C ratio vs. progenitor mass 
in the total ejecta of LIMS compared to observations in PNe 
(Palla et al. 2000). Stellar models are either from van den 
Hoek & Groenewegen (1997) (upper panel) or Ventura et 



T The values for three of the stars have been revised upwards 
later on by Onhaka et al. (2000) and found to be closer to those 
estimated by Lambert et al. (1986). 



Figure 3. In this figure we show the effect of lowering the extent 
of HBB in intermediate-mass stars on the predicted 12 C/ 13 C ra- 
tio. The dashed line shows the results of a model that adopts the 
standard yields of van den Hoek & Groenewegen (1997) (Model 
1); the continuous line shows the results of a model that accounts 
for the minimum HBB effect still compatible with the observa- 
tions. 

al. (2002) (lower panel), and have been computed for dif- 
ferent values of the initial stellar metallicity. Ventura et al. 
(2002) predict sistematically lower ratios than van den Hoek 
& Groenewegen (1997) in the mass range 3-6 Mq. However, 
their computations do not extend to the 1-3 Mq range, 
where model predictions could be usefully compared to ob- 
servations. This is the reason why we do not adopt the car- 
bon yields of Ventura et al. (2002) . Notice that in this figure 
we are comparing theoretical total yields with 12 C/ 13 C ratios 
observed in PNe. It would be better to compare the obser- 
vations with the yields at the tip of the AGB. We choose the 
total yields rather than those at the end of the AGB phase 
in order to compare the results from van den Hoek & Groe- 
newegen (1997) to those from Ventura et al. (2002). In fact, 
these latter authors do not give the yields at the tip of the 
AGB, but only the total ones. Theoretical 12 C/ 13 C ratios 
at the tip of the AGB from van den Hoek & Groenewegen 
(1997) are in general higher than those in the total ejecta. 
It is immediately seen that the adopted yields overestimate 
the 12 C/ 13 C ratio in the ejecta of low-mass stars. However, 
as demonstrated by Palla et al. (2000), as far as 12 C/ 13 C is 
concerned, the fraction of low-mass stars experiencing deep 
mixing does not affect significantly the overall results of the 
chemical evolution models, due to the fact that the Galac- 
tic evolution of 12 C and 13 C is mainly governed by stars in 
which this process is not expected to occur. Therefore, the 
main conclusions reached in the present paper should not 
be largely affected by the adoption of stellar yields taking 
extra- mixing into account. 

On the other hand, the discrepancy between model pre- 
dictions and observations can be solved if the adopted stel- 
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Figure 4. a) Temporal evolution of the oxygen isotopic ratio in the solar neighbourhood in the case of 17 production from only single 
stars (Model 1, dashed line), only novae (Model 2, dotted line) or novae plus single stars (Model 3, continuous line). The model le O/ 17 
ratios 4.5 Gyr ago are normalised to their solar values. The solar abundances by mass of ie O and 17 predicted by the models are 
also given at the bottom of the figure. The big dot is the observed meteoritic ratio; the vertical bar is the ratio observed in the local 
ISM (see text for references), b) Theoretical present-day oxygen isotopic ratio across the Galactic disc in the case of 17 production 
from only single stars (Model 1, dashed line), only novae (Model 2, dotted line) or novae plus single stars (Model 3, continuous line) vs. 
observations. Data are from Prantzos et al. (1996). 



lar yields overestimate the effect of HBB: running a model 
which adopts AGB yields computed with less HBB (tables 
from 22 to 31 of van den Hoek & Groenewegen 1997) leads to 
a 12 C/ 13 C ratio which decreases in the last 12 Gyr (Fig. 3), 
in agreement with observations. However, in this case the 
predicted solar abundance of 13 C is lower ( 13 C0 = 1.48 x 
10 -5 ). Notice also that the level of HBB required by this 
model is the lowest one still compatible with observations in 
stars. 

In Section 4.3 we examine in detail the possibility that 
novae contribute to the 13 C production. 



4.2 ir O production from single stars 

17 is produced by the NO cycle in both intermediate- and 
high-mass stars; contrary to 13 C, it is not produced in the 
1-2 Mq mass range where the NO cycle does not operate. 
It is mostly produced as a secondary element, owing to the 
dependence of its abundance on the initial 16 O content of 
the star. In Figs. 4a,b we show the predictions from Model 1 
( lr O production from single stars alone) as dashed lines. The 
nucleosynthesis prescriptions for 17 from LIMS and mas- 
sive stars are from Ventura et al. (2002) and Nomoto et al. 
(1997), respectively (see Table 1). It is worth noticing that, 
since the study of Nomoto et al. (1997) concerns only the 
He-cores and does not treat properly the products of the H- 
burning surviving outside the receding He-core, the yields of 
lr O from massive stars adopted by Model 1 are underesti- 
mated. Therefore, in this model 17 comes practically solely 
from intermediate-mass stars in the mass range 2.5-6 Mq. 



Nevertheless, in this model we have to reduce the O yields 
by Ventura et al. (2002) by a factor of ~ 6 in order not to 
overestimate the 17 solar abundance. In fact, Ventura et 
al. (2002) predict a too large production factor for lr O. This 
problem was already encountered by Prantzos et al. (1996), 
who adopted the 17 yields from LIMS by Marigo et al. 
(1996). They showed that the Marigo et al. yields have to 
be reduced by a factor of ~ 3 in order not to overproduce the 
solar 17 (cfr. figure If of Prantzos et al. 1996 and figure 7 
of Ventura et al. 2002). 

Model 1 predicts that the 16 0/ 17 ratio decreases dur- 
ing the first ~ 7 Gyrs of Galactic evolution, then it flat- 
tens out and finally increases at around t — 10 Gyr. At 
first glance, this is unexpected, since the 16 0/ 17 ratio is a 
typical primary to secondary ratio and hence, according to 
common wisdom, it should always decrease in time. The ex- 
planation of this behavior stems from the stellar mass range 
from which 17 originates in our model, and from the star 
formation history. The yields of Ventura et al. (2002) we 
adopt here are computed for stars in the mass range 2.5- 
6 M Q . The lifetime of a 6 M star is ~ 0.065 Gyr, while 
that of a 2.5 Mq star is ~ 0.75 Gyr. Therefore, stars in 
this mass range start restoring their newly synthesized 17 O 
rather early during Galaxy's evolution. An age of 7 Gyr (the 
age at which the oxygen isotopic ratio flattens out) corre- 
sponds to the lifetime of a 1.4 Mq star, while an age of 10 
Gyr (the age at which the oxygen isotopic ratio starts in- 
creasing) corresponds to the lifetime of a 1 Mq star. Stars 
with m < 2.5 Mq do not produce lr O in our model. From 
this fact and from an inspection of Fig. 5, it can be im- 



© 2003 RAS, MNRAS 000, 1-17 



10 D. Romano and F. Matteucci 




5 10 15 

Time (Gyr) 



Figure 5. The predicted temporal evolution of the total metal content (Z) and of the 16 0/ 17 ratio in the ISM for three different 
Galactocentric distances are shown as a thick solid and a thick dashed line, respectively. The star formation rate as a function of time 
is also shown, as a thin solid line. Its oscillating behavior at large Galactocentric distances is due to the presence of a threshold in the 
gas density below which the star formation stops (see Chiappini et al. 1997). 



mediately understood why, although 1 O is produced as a 
secondary element, the le O/ lr O ratio first flattens and then 
increases in time: this happens because a significant fraction 
of the stars which die at those times are not 17 O produc- 
ers. Moreover, they formed out of gas characterized by high 
16 0/ 17 ratios, and so the fraction of unburnt material that 
they restore into the ISM is characterized by the same high 
ratios. The predicted gradient across the Galactic disc can 
be easily explained by the same arguments. In particular, 
the high 16 0/ lr O ratio we find at inner radii at the present 
time is due to the formation of an important fraction of low- 



mass stars at early times, out of gas with a high le O/ 17 
content, which die at late times and restore into the ISM 
matter characterized by this high 16 0/ ir O ratio - and low 
metal abundance (Fig. 5, top panel). 

Summarizing: stars in the mass range 2.5-6 Mq are 
able to reproduce the solar abundance of 17 without re- 
quiring further ir O producers, but their theoretical yields 
have to be reduced (e.g., by a factor of ~ 6 in the case of 
the 17 O yields adopted in this work). However, once the so- 
lar lr O abundance is explained, theoretical predictions on 
both the temporal evolution of the 16 0/ 17 ratio in the 
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Table 3. Predicted and observed gradients (dex kpc 1 ) at the 
present time (*Gal = 14 Gyr). 



A( 12 C/ 13 C)/AR G A( 14 N/ 15 N)/AR G A( 16 0/ 17 0)/AR G 
Model 1 

4.6 79 
Model 2 

2.6 8 104 
Model 2s 

12 26 255 

Model 3 

2.7 8 108 
Model 3s 

4.3 26 226 

Model 3n 

4.8 19 164 
Model 3 m 

4.8 18 193 

Observed 

4.5 ± 2.2 a 19.7 ± 8.9 b 

a Boogert et al. (2000). 
b Dahmen ct al. (1995). 

solar neighbourhood and the 16 0/ 17 gradient along the 
Galactic disc are not in agreement with the observations. 

It is interesting to note that we reach the same con- 
clusion also if we adopt the metallicity dependent yields by 
Woosley & Weaver (1995) for massive stars and set them 
as the only lr O producers. In fact, by using the metallicity 
dependent yields by Woosley & Weaver (1995) for massive 
stars, the solar abundance of ir O is reproduced (within the 
errors), but the oxygen isotopic ratio remains almost flat 
in the past 4.5 Gyr. This result disagrees with previous in- 
vestigations suggesting that massive stars alone can explain 
the observed evolution of 16 0/ ir O in the Galaxy, once the 
metallicity effects on the stellar yields are properly taken 
into account (e.g., Timmes et al. 1995; Goswami & Prant- 
zos 2000). We argue that the differences could be due to the 
fact that in our models many stars in the mass range 1-2 
Mq die and restore matter with high 16 0/ 17 at late times. 
Moreover, the effect of the threshold in the star formation 
process can play a role in the same direction. 

In conclusion, our results suggest: i) the need for a ma- 
jor revision of the theoretical lr O yields from intermediate- 
mass stars on the one hand and ii) the occurrence of late 
lr O pollution from some different, long-lived stellar source 
on the other. 



4.3 CNO production during nova outbursts 

In previous papers (Romano et al. 1999, 2001) we adopted 
the 7 Li yields from novae by Jose & Hernanz (1998) and 
showed that novae are among the best candidates in order 
to explain the sudden increase of the lithium abundance in 
the ISM at [Fe/H] > —0.5 dex, as traced by dwarf stars 



in the solar vicinity hotter than 5700 K in the log e( 7 Li) - 
[Fe/H] diagram. The key ingredients which led to this result 
were i) the high overproduction factors of 7 Li with respect to 
its solar abundance and ii) the long delay with which novae 
restore their nuclearly processed material into the ISM. 

Large overproduction factors have been found also for 
13 C, 15 N and 17 (see figures 1 to 3 of Jose & Hernanz 
1998). These overproduction factors are so large that one 
can argue that novae can eventually account for most of the 
Galactic 13 C, 15 N and ir O (see the discussion in Section 1). 
Moreover, it is apparent that novae, restoring their newly 
synthesized 13 C, 15 N and 17 to the ISM on very long time- 
scales, act in the right direction in solving any discrepancy 
between model predictions and observations when dealing 
with the Galactic evolution of the CNO isotopic ratios. 

In Figs. 1, 4 and 6 we show as dotted lines the results 
of a detailed chemical evolution model in which all of the 
Galactic 13 C, 15 N and 17 are supposed to come from nova 
outbursts (Model 2, Table 1). In this model, the 13 C, 15 N 
and 17 production from novae is treated as primary, i.e., 
13 C, 15 N and 17 are assumed to be produced independently 
of the initial metallicity in the nova outburst. The average 
masses ejected in the form of 12 C, 13 C, 14 N, 15 N, ie O and 
17 assumed by the model are listed in column 5 of Table 1. 
We assume that the mean mass ejected during a single burst 
is ~ 2 x 10 -5 Mq and that each nova experiences 10 4 bursts 
during its life (see Section 3.2). Nucleosynthesis prescrip- 
tions for LIMS and massive stars are from van den Hoek & 
Groenewegen (1997 - their standard case in which the mass 
loss parameter for stars on the AGB, jjagb, is equal to 4 and 
does not vary with metallicity) and Nomoto et al. (1997), 
respectively. In this model, the 13 C, 15 N and lr O production 
from LIMS and massive stars is set to zero. It is found that 
it is necessary to lower the average nova yields of 13 C, 15 N 
and 17 from Jose & Hernanz (1998) in order to avoid to 
overestimate their meteoritic abundances. 

In Figs. 1, 4 and 6, panels a, the temporal evolution 
of the CNO isotopic ratios in the solar neighbourhood is 
displayed. If novae are the only sources of 13 C, 15 N and 
17 and produce them as primary elements, once the solar 
values are reproduced, the present-day local isotopic ratios 
fall to values smaller than observed (unless the 2 a level of 
uncertainty for measurements in the local ISM is assumed), 
due to the delayed contribution of too large amounts of 13 C, 
15 N and 17 from nova systems at late times. Moreover, 
the slope of the theoretical gradients is sistematically flatter 
than observed (Table 3) . 

If instead 13 C, 15 N and 17 from novae are secondary 
(Model 2s, Table 1), we obtain that: i) the decrease of the 
CNO isotopic ratios in the solar vicinity in the last 4.5 Gyr 
is more pronounced (Fig. 7, left panels); ii) the gradients 
steepen out (Fig. 7, right panels); in particular, in the case 
of nitrogen the agreement with the observations improves, 
while in the case of carbon it worsens (Table 3). 

It is interesting to analyse what is the effect of changing 
the free parameters involved in the computation of the nova 
outburst rate on the predicted CNO ratios. Changing the 
value of At, i.e., the delay time required to cool the WD at 
a level that guarantees strong enough nova bursts, from 2 
to 5 Gyr, results in a lower nova rate at the present time, 
and hence lower solar abundances for the CNO group ele- 
ments (both in the case in which the production of the minor 
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Figure 6. a) Temporal evolution of the nitrogen isotopic ratio in the solar neighbourhood. 15 N is only produced during nova outbursts. 
The model 14 N/ 15 N ratio 4.5 Gyr ago is normalised to its solar value. The solar abundances by mass of 14 N and 15 N predicted by the 
model are given at the bottom of the figure. The big dot is the observed meteoritic ratio; the vertical bars are the ratios observed in local 
diffuse (lower value) or dense (higher value) clouds (see text for references), b) Theoretical present-day nitrogen isotopic ratio across the 
Galactic disc. The dashed line are the actual model predictions (dotted line) offset to better compare the slope of the gradient with the 
data. 



CNO isotopes from novae is treated as primary and in the 
case in which it is treated as secondary), unless the value of 
an is properly increased")" . Increasing a so that the value of 
^novae(iGai) predicted with At — 5 Gyr is the same as that 
found with At — 2 Gyr results in flattening too much the 
gradients along the disc. 

It is clear that the nova contribution to the chemi- 
cal evolution of the Galaxy is at present only poorly con- 
strained: the rate of nova outbursts in the past is virtually 
unknown; the delay time At could vary from system to sys- 
tem; the mass ejected by a single nova during each burst 
is very uncertain (see D'Antona & Matteucci 1991 and the 
discussion in Section 3.2). However, our analysis allows us to 
remark at least two important things: i) first, if the fraction 
of WDs which enter the formation of nova systems is con- 
stant in time (similarly to what is assumed for the fraction of 
WDs entering the formation of binary systems ending up as 
Type la SNe), the nova yields of 13 C, 15 N and 17 computed 
by Jose & Hernanz for a solar composition of the matter ac- 
creted by the WD overestimate the solar abundances of 15 N 
and 17 0; ii) second, if novae are the only sources of 13 C, 
15 N and 17 0, the predicted 12 C/ 13 C, 14 N/ 15 N and 16 0/ 17 
ratios tend to decrease too fast from the time of Sun forma- 
tion up to now, especially in the case in which their yields 
are scaled with the initial stellar metallicity. 

In Figs. 1 and 4, we also show as continuous lines the 
predictions of a model (Model 3, Table 1) in which 13 C and 
17 are produced both by single stars and novae (as primary 

•f Notice that we cannot put constraints upon a and n singularly, 
but rather on the product an. 



elements, in this latter case). The decrease of the 12 C/ 13 C 
ratio in the solar neighbourhood in the last 4.5 Gyr can be 
now reproduced without requiring a reduction of the HBB 
strength in the van den Hoek & Groenewegen (1997) stan- 
dard set of yields. Also the gradient of this same ratio across 
the disc is well reproduced. However, in order not to over- 
produce the 17 abundance in the Sun, we need to lower in 
Model 3 the 17 O yields from both novae and intermediate- 
mass stars. In particular, in this model novae produce almost 
50 per cent of the solar 17 O, while the remaining 50 per cent 
comes from intermediate-mass stars. Actually, favouring a 
higher nova contribution with respect to that from LIMS 
would improve the agreement with the observations. How- 
ever, any percentage would be very tentative, due to the 
uncertainties in the nova and intermediate-mass star yields 
just mentioned before. Therefore, any acceptable solution 
would be by no means unique. 



4.4 Changing the nucleosynthesis prescriptions 
for massive stars and LIMS 

Finally, one could ask how firm are the results we presented 
in the previous section with respect to changes in the nucle- 
osynthesis prescriptions for novae, massive stars and LIMS; 
in other words, how much do our conclusions depend on the 
choice of the particular set of stellar yields we have adopted? 
In this section we show the results of a chemical evolution 
model, Model 3s, similar to Model 3, but where the yields 
from novae are scaled according to the initial metallicity of 
the star (Fig. 8, thin lines). In the same figure, we also show 
results from a model, Model 3n, similar to Model 3s, except 
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Figure 7. The theoretical behavior of the CNO isotopic ratios in the solar neighbourhood as a function of time (left panels) and across 
the Galactic disc at the present time (right panels) is shown for Model 2s as a dotted line compared to observations (see text and Figs. 1, 
4 and 6 for references). Model 2s is the same as Model 2, i.e., 13 C, 15 N and 17 are assumed to originate only from novae, but here the 
13 C, 15 N and 17 production from novae is assumed to be secondary rather than primary (see Table 1 and text for more details). The 
decrease of the CNO isotopic ratios in the solar neighbourhood in the last 4.5 Gyr is inconsistent with the observations, also if the 2 a 
level of uncertainty for measurements in the local ISM is assumed. 



for the fact that the yields for LIMS from van den Hoek 
& Groenewegen (1997) are computed with ?7agb increas- 
ing with increasing metallicity (from jjagb = 1 to jjagb = 
4). Furthermore, Model 3n adopts the yields from massive 
stars by Woosley & Weaver (1995) rather than Nomoto et 
al. (1997). It is worth noticing that the yields from LIMS 
increase with decreasing values of ??agb (i.e., smaller mass- 
loss rates) as a lower value of tjagb results in longer AGB 
lifetimes and therefore more thermal pulses (assuming that 



the amount of dredged-up matter during a thermal pulse 
is roughly constant). Moreover, the yields of 13 C and 17 
from Woosley & Weaver (1995) are larger than those from 
Nomoto et al. (1997), owing to the fact that they take into 
account nuclear burning outside the receding He-core. They 
also predict larger yields for 12 C, but slightly lower ones for 
16 0. 

As a consequence, Model 3n predicts a solar abundance 
of 12 C larger than that predicted by Model 3s. On the con- 
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Figure 8. The theoretical behavior of the CNO isotopic ratios in the solar neighbourhood as a function of time (left panels) and across 
the Galactic disc at the present time (right panels) is shown for Models 3s and 3n as a thin solid and a thick solid line, respectively, and 
compared to observations (see text and Figs. 1, 4 and 6 for references). As far as nova nucleosynthesis is concerned, Models 3s and 3n 
differ from Model 3 in the fact that the yields from novae are treated as secondary, i.e., they are scaled according to the initial metallicity 
of the star. Model 3n is the same as Model 3s, except for the nucleosynthesis prescriptions for single low-, intermediate- and high-mass 
stars (see Table 1). 



trary, the predicted solar abundance of 16 O is lower (see 
Table 2). Therefore, the secondary production of 13 C and 
15 N from novae increases, as the mean masses ejected from 
novae in the form of 13 C and 15 N are now scaled with a 
higher initial abundance of 12 C at each time. Owing to that, 
higher 13 C and 15 N solar abundances are found (Table 2). 
On the other hand, the secondary production of 17 from 
novae in the case of Model 3n is lower than in the case of 
Model 3s, since it scales with the initial abundance of 16 O 



in the star, which is now lower. However, Model 3n predicts 
for 17 a solar abundance higher than Model 3s, due to the 
non-negligible lr O contribution from massive stars when the 
yields of Woosley & Weaver (1995) are adopted rather than 
those of Nomoto et al. (1997). 

Once the CNO isotopic ratios are normalized to their 
theoretical solar values, it is seen that the model predictions 
for the 12 C/ 13 C and 14 N/ 15 N ratios from Model 3n, as far 
as both the temporal evolution at the solar ring and the 
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Figure 9. a) Temporal evolution of the oxygen isotopic ratio in the solar neighbourhood for two models in which 17 originates from 
intermediate- and high-mass stars as well as from novae (Models 3n and 3m, Table 1). The two models differ in the treatment of pristine 
17 O in stars: Model 3n (continuous line) assumes that all the 17 initially present in the gas out of which the star forms is destroyed 
later on by proton captures in the stellar interior, while Model 3m (long -dashed- short- dashed line) assumes that all the pristine 17 is 
preserved and ejected into the ISM at the death of the star, b) Theoretical oxygen isotopic ratio across the Galactic disc vs. observations 
for the same models as in panel a. 



present-day distribution across the disc are concerned, do 
not change appreciably with respect to those from Model 3s 
(Fig. 8 and Table 3). On the other hand, in the case of 
Model 3n the ie O/ lr O ratio shows a smoother decrease from 
the time of the birth of the Sun up to now, and a shallower 
gradient is predicted across the disc. This is due to the in- 
jection into the ISM of a non-negligible amount of 17 from 
massive stars. 

At this point, it has to be noticed that if the theoretical 
isotopic ratios are not normalized to the solar values pre- 
dicted by the models themselves, the situation is even worse. 
Our models predict solar abundances of 12 C, 14 N and 16 O in 
good agreement with the more recent abundance determina- 
tions in the solar photosphere by Holweger (2001) (except for 
12 C) and Allende Prieto, Lambert & Asplund (2001, 2002), 
but too low when compared to the meteoritic abundances by 
Anders & Grevesse (1989) (Table 2; see also Chiappini et al. 
2003). Since the isotopic ratios given in Anders & Grevesse 
can still be considered as the right Solar System values (P. 
Hoppe, private communication), it is likely that the solar 
abundances of the minor isotopes have to be scaled accord- 
ing to the revised values of the corresponding main isotopes. 
In other words, we should require an even lower contribu- 
tion from novae and/or LIMS to the Galactic 13 C and 17 
abundances than suggested above. This may be a precious 
information for people dealing with stellar evolution and nu- 
cleosynthesis. 

Finally, we address the problem of 17 destruction in 
stars. 17 is destroyed in stellar interiors by proton captures 
through the reactions 17 (p, a) 14 N and 17 (p, 7) 18 F (this 
latter occurs only at high temperatures). All the models 



previously discussed assume that all the O injected into 
the ISM by stars of different masses is newly produced, i.e., 
all the lr O present in the protostellar nebula is destroyed 
in the hot stellar interior. In Fig. 9a,b we compare results 
from Model 3n (solid lines) to those from Model 3m (dashed 
lines). Model 3m is the same as Model 3n, except for the 
fact that now all the pristine lr O is assumed to survive and 
to be returned back into the ISM at the death of the star. A 
realistic situation should be probably something in between 
the two. Assuming that all the 17 O present in the gas out 
of which the stars form is preserved, causes a steepening of 
the 16 0/ 17 gradient and a more pronounced decrease of 
the ratio in the solar neighbourhood in the last 4.5 Gyr. 
Obviously, in the case of Model 3m we need also to further 
reduce the production of newly formed 17 O in order not to 
overestimate its solar abundance. We choose to lower the 
yields of 17 from intermediate-mass stars (see Table 1). 



5 CONCLUSIONS 

In this paper we discuss the problem of the evolution of 
the CNO isotopes in both the solar vicinity and the disc of 
the Galaxy. In particular, we analyse the implications of up- 
dated results from stellar nucleosynthesis studies (including 
nova nucleosynthesis). Our main conclusions can be summa- 
rized as follows: 

(i) A model which con- 

siders only single low-, intermediate- and high-mass stars 
as 13 C, 15 N and ir O producers, is in trouble with repro- 
ducing the data relevant to the temporal variation of the 
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CNO isotopic ratios in the solar neighbourhood as well as 
their behavior across the Galactic disc at the present time, 
unless one does not allow for a proper revision of current 
stellar yields. In particular, in the specific case of the car- 
bon isotopic ratio, we show that a good agreement with the 
observations can be obtained only by reducing the strength 
of HBB in intermediate-mass stars. 

(ii) Current theoretical yields of lr O from intermediate- 
mass stars are overestimated; this may suggest the need of 
considering some neglected reactions resulting in 17 O de- 
struction inside stars. 

(iii) It has been claimed that novae can contribute to large 
amounts of 13 C, 15 N and ir O and eventually explain the 
whole solar abundances of these elements (e.g., Starrfield et 
al. 1972, 1974; D'Antona & Matteucci 1991; Matteucci & 
D'Antona 1991; Woosley et al. 1997; Jose & Hernanz 1998). 
In this paper, we adopt detailed nucleosynthesis in the ejecta 
of classical novae as published by Jose & Hernanz (1998) for 
a grid of hydrodynamical nova models spanning a wide range 
of CO and ONe WD masses (0.8 - 1.35 M ) and mixing lev- 
els between the accreted envelope and the outermost shells 
of the underlying WD core (25% -75%). We find that, when 
included in a detailed model for the chemical evolution of the 
Milky Way, they produce 12 C/ 13 C, 14 N/ 15 N and 16 0/ 17 
ratios decreasing with increasing metallicity, i.e., decreasing 
with time at the solar radius and increasing with Galac- 
tocentric distance at the present time, in agreement with 
the trends inferred from observations. However, if novae are 
the only 13 C, 15 N and 17 producers, the CNO isotopic ra- 
tios are found to decrease too steeply at ii = Rq, often 
producing present-day ratios across the Galactic disc lying 
at the lower observational limit suggested by the observa- 
tions at each radius. These results are almost independent 
of whether the production of 13 C, 15 N and ir O from novae is 
treated as primary or secondary. In fact, also in the case of 
primary production, 13 C, 15 N and ir O behave as secondary 
elements from a point of view of Galactic chemical evolu- 
tion, owing to the long delay with which they are restored 
into the ISM. However, in the case of purely primary pro- 
duction (from the nucleosynthesis point of view) of 13 C, 15 N 
and ir O from novae, natter gradients and less pronounced 
decreasing trends at the solar radius are found for all the 
ratios. 

(iv) A model in which 13 C and 17 O are produced by 
intermediate- and high-mass stars as well as novae fits better 
the observations; however, in the case of 17 O it is necessary 
to lower by hand the yields from both intermediate-mass 
stars and novae in order to reproduce the solar lr O abun- 
dance. We know that 17 O is always produced as a secondary 
element in intermediate- and high-mass stars. However, our 
models do not rule out that at least a minor fraction of the 
solar 17 might be of primary origin. This fraction would 
originate from nova systems. 

(v) The behavior of the 14 N/ 15 N ratio along the Galac- 
tic disc and as a function of time in the solar neighbour- 
hood seems to suggest that 15 N has been produced on long 
time-scales, even longer than those of 14 N. Novae are the 
best candidates for producing 15 N on very long time-scales 
in the framework of the presently available nucleosynthesis 
calculations. However, large uncertainties are still present in 
the 15 N abundance determinations in the local ISM and this 



prevents us from drawing firm conclusions on the origin of 
this element. 
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